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Abstract 

The connection between the method of comparison equations (generalized WKB 
method) and the Ermakov-Pinney equation is established. A perturbative scheme 
of solution of the generalized Ermakov-Pinney equation is developed and is applied 
to the construction of perturbative series for second-order differential equations with 
and without turning points. 
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1 Introduction: the general case 



The WKB method is one of the main tools in quantum mechanics [1,2] and 
quantum field theory [3]. From the mathematical point of view this method 
coincides with the construction of the asymptotic expansions for the solution to 
the second-order differential equations with the small coefficient of the second 
derivative term [4,5]. In this note we consider an interesting modification of 
the WKB method - the so called method of comparison Eqs. [6,7] and analyse 
its relation with the Ermakov - Pinney equation [8,9]. Let us consider the 
second-order differential equation 
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^) = y {x)u{x) , (1) 

where e is a small parameter. Let us suppose that one knows the solution of 
another second-order differential equation 



-JLl = -Sl\c)U(a) . (2) 
In this case one can represent an exact solution of the Eq. (1) in the form 



where the relation between the variables x and a is given by the equation 



On solving Eq. (4) for a(x) and substituting it into Eq. (3) one finds the 
solution to Eq. (1). Eq. (4) can be solved by using some iterative scheme. 
The application of such a scheme is equivalent to the application of the WKB 
method or, in other words, to the construction of the uniform asymptotic ex- 
pansion for the solution of Eq. (1). The method of construction of the solution 
to Eq. (1) by means of the solutions of Eqs. (2) and (4) is called the method of 
comparison equations and the function Q 2 (cr) is called the comparison func- 
tion [7]. 

The iterative scheme of solution of Eq. (4) depends essentially on the form of 
the comparison function fi(cr). A reasonable approach consists in the elimina- 
tion of a from Eq. (4) and its reduction to a form, where the only unknown 
function is (da/dx). The equation obtained for the case of the simple com- 
parison function Q 2 (a) = 1 coincides with the Ermakov-Pinney equation, this 
will be explicitly shown in the next section. For more complicated forms of 
the comparison function Q 2 (a) useful for the description of various physical 
problems starting with the classical problem of tunneling for non-relativistic 
quantum mechanics [1,2,10] and ending with the WKB-type analysis of the 
linearised equations for cosmological perturbations [11,12], the corresponding 
equation will acquire a more intricate form and we shall call it a generalised 
Ermakov-Pinney equation. The role of an unknown function in this equation 
is played by the function 



-1/2 

y(*) = ( ^ 1 , (s) 
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while the equation itself has the form 



v ( dy\ ( dy d 2 y d 3 y\ 

where the concrete forms of the functions K and L will be defined in the next 
sections. We shall represent the function y(x) as a series 



y = E Vnt n , (7) 

n=0 

and the general recurrence relation connecting the different coefficient func- 
tions y n can be represented as 



'd n K d n ~ x V 



de n de n - 1 



0. (8) 



e=0 

To implement this formula one can use the combinatorial relation [13] 



= (« + «)<">, (9) 
dx n K ' w 

where on the right-hand side one has a binomial expression wherein the powers 
are replaced by the derivatives. When the function y(x) is constructed up to 
some level of approximation one can, in principle, find a(x) and substitute for 
it into Eq. (3). 

We shall consider the implementation of this scheme for four choices of the 
comparison function: f2 2 (<r) = l,Q 2 (a) = a,Q 2 (a) = exp(cr) — 1 and Q 2 (a) = 
a 2 — a 2 . At the end of the paper we shall discuss the relation between the 
different versions of the WKB method and equations of the Ermakov-Pinney 
type. 



2 Differential equation without turning points 



In the case for which the function uu 2 (x) does not have zeros (turning points) 
it is convenient to choose the comparison function Q 2 (a) — 1. In this case the 
function U(a) is simply an exponent, while Eq. (4) can be rewritten as 



d 2 y(x) 9/ , . . 1 . , 

^ = " <^> - m • < 10 > 
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where the function y(x) is defined by Eq. (5). The above equation is nothing 
but the well known Pinney or Ermakov-Pinney Eq. [8,9], which can be solved 
perturbatively with respect to the parameter e [14]. We shall give here some 
general formulae for such a solution. It is convenient to rewrite Eq. (10) as 



Luy-l=ey 3 y, (11) 

where a "dot" denotes the derivative with respect to x. We shall search for 
the solution of Eq. (11) in the form (7). The zero-order solution of Eq. (11) is 



Vo = ^~ 1/2 , (12) 
and the general recurrence relation following from Eq. (11) for n > 1 is 



y n—1 n—1 n—1 n—1 



y n = - -7-3 J2 ^ 5[n- (k 1 + k 2 + h + k 4 )\ yk 1 yk 2 yk 3 yk4 

4 2/0 fci=0fc 2 =0fc 3 =0fc 4 =0 



■y n—1 n — 1 



+ T^n, Z ■ ■ ■ ^ 5[n - l - + k 2 + h + k 4 )] y kl y k2 yk 3 yk 4 , (13) 
4u; yo kl =o k 4 =o 

where S() is the Kronecker delta symbol. It is now easy to obtain from the 
general expression (13) expressions for particular values of y n 



%i . y\ . 32/o2/i nAU , 

y^ = ~~A — h + a 9 ( 14b ) 

4y 4u 2 4u 2 y 
yf 3J/22/1 2/2 

S/3 = --2 - — — + 



2/0 yo 4cj 2 

32/Q2/2 32/i2/i 32/o2/f 
4cj 2 2/o 4w 2 2/o 4cu 2 2/ 2 ' 

Let us end by noting that the standard WKB approximation corresponds to 
the case discussed in this section and is of course only valid away from turning 
points. 



3 Turning point: the Langer solution 



In the case of the presence of a linear zero in the function u 2 (x) one can 
use the Langer solution [10]. In terms of the method of comparison equations 
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it means that one chooses the comparison function Q 2 (cr) = a. In this case 
Eq. (4) becomes 



2 , , (daV ida\ 1 ' 2 d? (day 11 ' 2 

Dividing this equation by (da/dx) 2 and differentiating the result with respect 
to x we get an equation which depends only on the derivative a and not on 
the function a. Such an equation can be rewritten in the form 



2uuy % + 4c^ 5 = 1 + e (3yyy 4 + y<V) , (16) 

where the function y(x), as in the preceding section, is defined by Eq. (5). 
Again we shall search for the solution of Eq. (16) in the form (7). The equation 
for y is 



2djujy 6 + AuYvl = (coYo) ' vl = 1 ■ (17) 
This equation my be rewritten as 



Wo y x (^o) = <" > ( 18 ) 
and integrated 



Finally 



Vo = - 



"-(f/^f. (19) 



1/3/- \s 

V2[2j udx ) ■ (20) 



Comparing the terms arising in Eq. (16) as coefficients of e n with n > 1 we 
have 



AuYoVn + (l2uuy 5 + 20u 2 y o yo) Vn = ^~ 1/2 (ooYoVn) ' = F n , (21) 

where F n contains the terms which depend on y k , < k < n — 1, and can be 
written as 
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n— 1 n—1 



F n = J2 • ' • { S [n-ih + h + h + h + h + h)} (^uj 2 y kl y k2 - 2uuy kl y k2 ) 

fci=0 k 6 =0 

+5[n-l-(k 1 +k 2 +k 3 +k 4 +k 5 +k 6 )\ (3y kl y k2 + yffy k2 ) y k3 y ki y k5 y k(i } ■ (22) 
The general expression for y n is 

y n = -^Jdxu>F n . (23) 
Let us also give the explicit expressions for the first terms of the expansion (7) 

1 t 

Vi = J dxuj (3y y y* + y^y^) (24a) 

1 r 

Vi = J dxuJ (-WuuylyQ - 20uj 2 y 1 y 1 y* - A0uj 2 y y 2 yQ 

+ zymvo + ^Mivo + y { ?vl + %o 3) 2/i2/ 4 ) (24b) 

1 t 

2/3 = 47j3^5 J dxuJ (-^^yiVo ~ ^Ouuyfyl - 20u 2 y 2 yiyo 
- 80u 2 y y 2 yiyl + 3y 2 y yo + Sy y 2 y^ + l2y Q y Q y 2 yl + Zyiy\vt 

+ ^ymyiVo + Hyoymy + ^mmyiyo + y\ yo 

+ 5y^y 2 y 4 + hy { ?mj$ + IQy^yho) . (24c) 

Before closing this section we may compare our results with some formulae, 
obtained by Dingle [7]. In the vicinity of the turning point, which we choose 
as x = 0, the function u 2 (x) can be represented as 



^)=7i* + f* 2 + |* 3 + |* 4 + --- (25) 
while the function a(x) looks like 



a(x) = a + viz + ^x 2 + ^x 3 + ^ + ■ ■ ■ (26) 

In reference [7] the first five coefficients <7o, a±, a 2 , a 3 and 04, as functions of 
the coefficients 71, 72, 73 and 74, were obtained in the lowest approximation by 
using a system of recurrence relations. We can reproduce all these formulae 
by only using the general formula (20). Indeed, on writing 



da _ I _ u(x) ^ 



dx 



,2 



1/3 
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using expansions (25) and (26) and comparing the coefficients of different 
powers of x we get the formulae 



O"0 


= 






o\ 


1/3 
= 7i 








1 _ 
= 5 71 


2/3 

72 






1 _ 
= 7 71 


2/3 

73 - 


12 
175 


(7 4 


1 _ 
= 9 71 


2/3 

74 - 


44 
315 



-5/3 2 

7i 7 2 

-5/3 148 _ 8 /3 3 , . 

7i 7273 + 7^71 7 2 , (28) 

which coincide with those obtained in [7] . The value of <7o comes directly from 
Eq. (15) ignoring the e-term and using the relations (25) and (26). 

We can also find the e-dependent corrections to the function a(x) in the neigh- 
bourhood of the point x = by using the results (28) and the recurrence re- 
lations (24a) etc. We shall only give here the first correction, proportional to 
e, to the coefficients a and U\. The value of a comes directly from Eq. (15) 
using the relations (25), (26) and (28) 



a = e 



7s7i 



5/3 



14 



97 2 2 7i 
140 



(29) 



In order to obtain o\ we introduce into the definition of the function y(x), 
Eq. (5), the expansion (26) and 



-1/2 



'3cr% 0-3 

2(7/" ' V8of 4(1! 



1 a2 

1 — - — x+ 



2 , / 30" 2 0"3 °"4 
X + 



5(7 U X* 



8(7i 2 12o"i 16a? J 



+ 



(30) 



Now, on using the recurrence formula (24a) we find the leading (constant) 
term in the expression for y 1 



1 / 1 3 -7/2 . 1 -5/2 C4C1 

yi = - -^1 +2<wi 



■3/2 % 



+ 



(31) 



From the formula 



da 



dx y 2 (y + eyif 



(32) 



it is easy to find the first correction to the first coefficient <7i in the expan- 
sion (26) 
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1/3 Hi 



(33) 



Substituting into Eq. (33) the leading term of y\ from Eq. (31), of yo from (30) 
and using the explicit expressions for the coefficients <7j from Eq. (28) we finally 
obtain 

1/3 . / 7 3 _ii/3 7 _8/3 . 1 -5/3 \ (OA s 

ai = 71 + 6 l 225 7271 " 135 737271 + 54 7471 J • (34) 

The results (29) and (34) was found by Dingle [7], but in his article there were 
two sign errors in the second line of equation (69). 



4 Turning point: more complicated comparison function 

Let us consider the differential Eq. (1), when the function uj 2 (x) has a linear 
turning point, but instead of Langer's comparison function Q 2 (a) = a consider 
a more complicated comparison function, also having linear turning point, 
namely 

=e a -1 . (35) 
Now, Eq. (4) will have the form 

Isolating the function e a , taking its logarithm and differentiating it with re- 
spect to x we arrive to the equation which involves only the derivative a. This 
equation can be written down as 

uY - 2uuy 6 - 4u 2 yy 5 + 1 = e (yy 3 - 3yyy 4 - y {3) y 5 ) . (37) 
To find the function y one should solve the equation 

coY - 2ojojy$ - AooYvl + 1=0. (38) 
On introducing the new variable 

z = ujYo (39) 
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one can rewrite Eq. (38) as 

. _ (z+ l)u 
It is also convenient to introduce the variable 

v = \fz . 
Eq. (40) now becomes 



2v 2 i) 

= UJ 



1 + v 2 

which can be integrated, obtaining 

2v — 2 arctan v = J cudx . 
Finally, for y we have the implicit representation 

2 1 2 1 f . 

y n arctan ooy n = — / uodx . 

uj 2,co J 

We may now write the equation defining the recurrence relation for ;, 
where 



n—1 n—1 

F n = -LU 2 E 5 N~( /i; l + /i; 2 + ^3 + ^4)]l/fc 1 l/fe 2 l/fc3l/fc4 

fe 1= fc 4 =0 
n—1 n—1 

+2cjw "' E 8[n-(h+fa+fa+fa+k 5 +k 6 )]y kl - ■ -y ke 

fci=0 fc 6 =0 
n—1 n—1 

+4w 2 E "' E S[n-(k 1 +k 2 +k 3 +k4+k 5 +kfi)]y kl y hll '-'y ko 

fci=0 fe 6 =0 
n—1 n—1 

+ E E d[n- 1 -(k 1 +k2+k 3 +k 4 )]y kl y k2 y k3 y k4 

kl =0 fc 4 =0 
n—1 n—1 

~ 3 E "' E 5 ["-!- ( fc i +^2+^3+^4+^5+^)] y^ykzVks ■■ -yk 6 

kl =0 k 6 =0 
n—1 n—1 

"E ••• E ^^-l-l^i + ^ + ^ + ^ + ^ + ^e)]^^---^ • (46) 

fci=0 fc 6 =0 

The solution for y n is 

l —JdxuGF n) (47) 



where 

«„=,(-/ 



(48) 



yl{x')) ' 

Here we give the explicit expressions for the first two functions F\ and F 2 

Fi = myl - 3Moi/o - yo ] yt , ( 49a ) 

F 2 = -5u 2 y 2 y 2 + SOuxvyfyo + 20^^^ + 40uj 2 y y 2 y% + y\y^ 

+ $ymyt - Zymyo - 3Mi2/o - l2 mymyl - yi 3) 2/o - %o 3) 2/12/0 • ( 49b ) 

5 Two turning points 



Let us now consider the differential Eq. (1), when the function u 2 (x) has two 
turning points. This situation, describing the tunneling through a potential 
barrier was considered in [6,2]. In this case the comparison function can be 
chosen as 



n 2 (a) =a 2 -a 2 . (50) 
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The solutions U(a) of Eq. (2) are the well-known parabolic functions [2]. 
Substituting Q 2 (a) of Eq. (50) into Eq. (4) one can isolate a and after the 
subsequent differentiation one can write down the equation for the function 

y(x) 



2^ V + a 2 + eyy 3 = 2uuy & + 4uj 2 yy 5 + e {Zyyy A + y% 5 ) . (51) 
The lowest-approximation y can be found from the equation 



2 Ju 2 y 4 + a 2 = 2djuy b + 4uj 2 yy 5 . (52) 



This equation can be integrated and the solution y can be written down in 
the implicit form 



uoyl arcsinh^^ = — [ udx . (53) 

uj a uj J 



To get the recurrence relations for the functions y n , n > 1 it is convenient to 
take the square of Eq. (51) 



4dj 2 uj 2 y 12 + lQuj 4 y 2 y w + 16djio 3 yy u - 4cjV - 4a 2 = 

e{4yy 3 - 12uuyyy 10 - 2Auj 2 yy 2 y 9 - Auujy {3) y u - 8uj 2 y {3) yy w ) 

-e 2 (9fy 2 y 8 + y^ 2 y 10 + 6y^yyy 9 ). (54) 



On comparing the terms containing the n th powers of the small parameter e 
one obtains 



fy n + 9Vn = K , (55) 
where 
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f = 32uj%y 1 ° + 16djuj 3 y 1 1 (56a) 
g = ASuWy] 1 + lQOu 4 y 2 y 9 + 176^ 3 y y 10 - 16uj 2 y 3 (56b) 

n— 1 n—l 

F «= E ••• E 5[n-(A; 1 + --- + A; 12 )] 

fci=0 fei2=0 

x 4 (uj 2 uj 2 y kl y k2 + lQuj 4 y kl y k2 + lQdjuj 3 y kl y k2 ^ y hi ■ --y kl2 

n—l n—l 

- Au 2 ■ • • E 5 l n ~ H h Ml • • • yk 4 

kl =0 fc 4 =0 

{n—l n—l 
~ 4 E ' • ' E 5 ( n ~ 1 ~ ( fc i + ' • ' + yk 1 yk 2 yk 3 yk 4 
kl =0 fc 4 =0 

n—l n—l 

+ E • • • E 5 [« - 1 - ( fc i + • • • + ^12)] 

fc 1= k 12 =0 

x (l2uuy kl y k2 y k3 + 24uj 2 y kl y k2 y k3 

+4uujy (3) y k2 y k:i + 8u 2 y { k ^y k2 y k3 ) y kA ■ ■ -y kl2 ] 

n-2 n-2 

-e 2 ]T ■■■ E 5[n-2-(h 1 + .'. + k 12 )] 

kl =0 fci2=0 

x (9y^y k2 y k3 y ki + y k ?y k2 y k3 y ki + Qy { k ^yk 2 yk 3 yk 4 ) j/k 5 • • • y kl2 ■ (56c) 
The integration of Eq. (55) gives 



Vn = h I dx ^~r ' ^ 



where 



G = exp ^ J dx'j^J . (58) 

The explicit expressions for the coefficients F n are rather cumbersome and we 
shall only write down the coefficient F 1 

F l = A^uWy^l 1 + 32^02/0° + W&>Hlvivl 

+ lQdjuj^y] 1 + mdj^yoy^ 1 ? - lQu 2 y iy 3 + 4y y 3 

- Udjcoyoyoyl ~ ~ 4wuy^Vo - 8u 2 y^y y] . (59) 
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6 Discussion: WKB method and the generalized Ermakov-Pinney 
equation 



In this paper we have studied the relation between the WKB method for 
the solution of second-order differential equations and the Ermakov-Pinney 
equation. For some non-trivial applications of the WKB method it is found 
that instead of the standard form of the Ermakov- Pinney equation one is lead 
to its generalizations which are rather different to the known generalizations 
of Ermakov systems. 

Let us note that although the Ermakov equation was already reproduced in 
the WKB context (first of all in the important paper by Milne [15]), its gen- 
eralizations arising in the process of application of the comparison function 
method [6,7] were not yet studied, at least to the best of our knowledge. In 
this concluding section of our paper we shall try to give a short review of 
the history of the Ermakov-Pinney equation and its applications to quantum 
mechanics and to similar WKB-type problems and discuss the physical sig- 
nificance of its non-trivial generalization (for a more detailed account of the 
history of the Eramkov- Pinney equation see e.g. [16]). 

The relation between the second-order linear differential equation 



dx 2 



+ u 2 (x)u(x) = 



(60) 



and the non-linear differential equation 



d 2 p 2 a 

d^ + uj{x)pix) = 7(x)' 



(61) 



where a is some constant was noticed by Ermakov [8], who showed that any 
two solutions p and u of the above equations are connected by the formula 



r P n, 

'Gi— - a . 
u z 



(62) 



The couple of Eqs. (60) and (61) constitute the so called Ermakov system. An 
important corollary was derive [8] from the formula (62). Namely, on having 
a particular solution p{x) of Eq. (61) one can construct the general solution 
of Eq. (60) which is given by 



u(x) = c±p(x) exp 



dx 
p 2 (x) 



+ c 2 p(x) exp 




(63) 
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It is easy to see that the function p(x) plays of the role of an "amplitude" of the 
function u(x), while the integral / [dx/p 2 (x)} represents some kind of "phase". 
Thus, it is not surprising that the Ermakov equation was re-discovered by 
Milne [15] in a quantum-mechanical context. On introducing the amplitude, 
obeying Eq. (61) and the phase, Milne constructed a formalism for the solution 
of the Schrodinger equation which was equivalent to the WKB method. Milne's 
version of the WKB technique was extensively used for the solution of quantum 
mechanical problems (see e.g. [17]). In the paper by Pinney [9] the general form 
of the solution of the Ermakov equation was presented, while the most general 
expression for this solution was written down by Lewis [14]. The equation, 
sometimes called Ermakov-Milne-Pinney equation, has also found application 
in cosmology [18]. 

A very simple physical example rendering trasparent the derivation of the 
Ermakov- Pinney equation and its solutions was given in the paper by Eliezer 
and Gray [19]. The motion of a two-dimensional oscillator with time-dependent 
frequency was considered. In this case the second-order linear differential 
Eq. (60) describes the projection of the motion of this two-dimensional oscil- 
lator on one of its Cartesian coordinates, while the Ermakov-Pinney Eq. (61) 
describes the evolution of the radial coordinate p. The parameter a is nothing 
more than the squared conserved angular momentum of the two-dimensional 
oscillator. Thus, the notion of the amplitude and phase acquire in this example 
a simple geoemtrical and physical meaning. 

An additional generalization of the notion of Ermakov system of equations 
was suggested in the paper by Ray and Reid [20]. Thy consider the system of 
two equations 



where g and / are arbitrary functions of their arguments. The standard Er- 
makov system (60), (61) corresponds to the choice of functions g(£) = and 
/(£) = a C One can show that the generalized Ermakov system (64), (65) has 
an invariant (zero total time derivative) 




(64) 



and 




(65) 



1 



\p J \uj \ ax ax 



) 



2 



2 



(66) 
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where 



(67a) 




(67b) 



The invariant (66) establishes the connection between the solutions of Eqs. (64) 
and (65) and sometimes allows one to find the solution of one of these equa- 
tions provided the solution of the other one is known (just as in the case of 
the standard Ermakov system). However, for the case of an arbitrary couple 
of functions / and g a simple physical or geometrical interpretation of the 
Ermakov system analogous to that given in [19] is not known. 

In the present paper we, on one hand, have established the connection be- 
tween the generalized WKB method (method of comparison Eqs. [6,7]) and 
the Ermakov - Pinney equation and on the other hand we have obtained a 
generalization of the Ermakov - Pinney equation, different to those studied in 
the literature. Let us briefly summarize our approach. One has two second- 
order differential Eqs. (1) and (2) and the solution of one of them (2) is well 
known and described. Using the previous convenient analogy [19], one can say 
that in this case one has two time-dependent oscillators, evolving with two 
different time-parameters. One can try to find the solution of Eq. (1) by rep- 
resenting it as a known solution of Eq. (2) multiplied by a correction factor. 
This correction factor plays the role of the prefactor in the standard WKB 
approach while the known solution of Eq. (2) represents some kind of general- 
ized phase term. Further, the prefactor is expressed in terms of the derivative 
between two variables x and o (or in terms of the oscillator analogy, between 
two times). On writing down the equation definining this factor, which we 
have denoted by y(x) (see Eq. (5)) we arrive to the Eq. (4) which could be 
called generalized Ermakov-Pinney equation. For the case when the function 
uj 2 (x) does not have turning points the comparison function Q 2 (a) can be 
chosen constant (the second oscillator has a time-independent frequency) and 
the the equation defining the prefactor becomes the standard Ermakov-Pinney 
Eq. (10). In terms of the two-dimensional oscillator analogy [19] this means 
that we exclude from the equation for the radial coordinate the dependence 
on the angle coordinate a by using its cyclicity, i.e. the conservation of the 
angular moment. 

In cases for which the function u 2 (x) has turning points, as in Sees. 3, 4 
and 5, the comparison functions are non-constant and instead of the standard 
Ermakov-Pinney equation we have its non-trivial generalizations. The common 
feature of these generalizations consists in the fact that the corresponding 
equations for the variable y depend explicitly also on the parameter a which 
is not excluded automatically. To get a differential equation for y, one should 
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isolate the parameter a and subsequently differentiate with respect to x. As 
a result one gets a differential equation of higher-order for the function y(x). 
Remarkably, for the perturbative solution of these equations one can again 
construct a reasonable iterative procedure. 

Again, it is interesting to look at the generalized Ermakov-Pinney equation as 
an equation describing a two-dimensional physical system in the spirit of the 
reference [19]. In our case one has 



d 2 y(x) 9/ . . . H 2 (a) 

-ir = 7^. < 68 > 

where y plays the role of a radial coordinate, a resembles a phase or an angle 
(e.g. the position of a hand of a clock) and x is a time parameter. It is impor- 
tant to notice that according to the definition of the variable y (5) there is a 
relation between the radial and angle coordinates 



da 1 . . 

& = (69) 

which on interpreting a as a phase corresponds to constant (unit) angular mo- 
mentum. In our case however the right-hand side of the radial Eq. (68) contains 
an explicit dependence on the "angle" a. Thus, the couple of Eqs. (68) and 
(69) cannot be represented as a system of equations of motion correspond- 
ing to some Lagrangian or Hamiltonian as in Ref. [19]. Perhaps, this system 
could be described in terms of non-Hamiltonian dissipative dynamics, but this 
question requires a further study. 
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